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Abstract. In a minimal supersymmetric extension of the standard model (MSSM) with explicit 

CP violation a mass of the lightest Higgs boson H\ in the range 7 GeV < Mh 1 ^ 10 GeV is 
experimentally allowed by present accelerator limits. In the same scenario a lightest neutralino as 
light as 2.9 GeV can be a viable dark matter candidate, provided that a departure from the usual 
GUT relation among gaugino masses is assumed. 

PACS. 12.60.Jv, Supersymmetric models - 14.80.Cp Non-standard-model Higgs bosons - 95.35.+d 
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1 Introduction 

In SUSY models, CP violating phases in the soft terms 
can considerably enrich the phenomenology without 
violating existing constraints. In particular, in such a 
scenario it is well known that the mass of the lightest 
Higgs boson H\ can be much lighter compared to the 
CP conserving case. We wish here to take a fresh look 
at how light rriff 1 can be, and discuss the cosmolog- 
ical lower bound to the mass of the relic neutralino 
X, when standard assumptions are made for the origin 
and evolution of its relic density. For more details of 
this analysis, see PQ. 

2 MSSM with explicit CP violation: the 
CPX scenario 

In the presence of sizable CP phases in the relevant soft 
SUSY breaking terms, a significant mixing between the 
scalar and pseudo-scalar neutral Higgs bosons can be 
generated in the Higgs potential at one loop [2]. As 
a consequence of this, the three neutral MSSM Higgs 
mass eigenstates, labeled in order of increasing mass 
as Mh 1 < Mh 2 < Mh 3 , have no longer definite CP 
parities, but become mixtures of CP-even and CP- 
odd states. Due to the large Yukawa couplings, the 
CP-violating mixing among the neutral Higgs bosons 
is dominated by the contribution of third-generation 
squarks and is proportional to the combination: 
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with f — t,b. Here fi is the Higgs-mixing parameter 
in the superpotential and A t denotes the trilinear soft 
coupling. 

In presence of CP violation, the mixing among neu- 
tral Higgs bosons is described by a 3 x 3 real orthogonal 
matrix O: 

(0 1 ,0 2 ,a) T = O(!/ 1 ,H 2 ,!/ 3 ) T , (2) 

where the elements O^i and O^i are the CP-even 
components of the i-th Higgs boson, while O a i is the 
corresponding CP-odd component. 

The Higgs-boson couplings to the SM and SUSY 
particles can be modified significantly due to the CP 
violating mixing. Among them, one of the most im- 
portant ones is the Higgs-boson coupling to a pair of 
vector bosons, gHiW i which is responsible for the pro- 
duction of Higgs bosons at e + e~ colliders: 

Chw = gM w (w+w-» + -hrZ^zA Y^g^wH, , 

V Zc w J i=1 

(3) 

where 

9HiW = cp Ofai + sp 0<j, 2 i , (4) 

when normalized to the SM value. Here we have used 
the following abbreviations: sp = sin/3, cp = cos/3. 
tp = tan/3, etc. We note that the two vector bosons W 
and Z couple only to the CP-even components 0^ 11 2 i 
of the i-th Higgs mass eigenstate, and the relevant cou- 
plings may be strongly suppressed when the i-th Higgs 
boson is mostly CP-odd, 2 al ~ 1 » 0% xl , 0%^. 

The so called CPX scenario is defined as a show- 
case benchmark point for studying CP- violating Higgs- 
mixing phenomena [3]. Its parameters are all defined 
at the electro-weak scale, and are chosen in order to 
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enhance the combination in Eq.JT]). In this scenario, 
SUSY soft parameters are fixed as follows: 



> 



M f 



M f 



M Q 3 ~ J "C/ 3 _ 1V± D 3 

M L 3 = m e 3 = A/susy; \h\ = 4A/ SUS Y, 
|A t , 6ir |=2Ms USY ,|M 3 | = l TeV, (5) 

where, with a usual notation, Q, L, U, D and E indi- 
cate chiral supermultiplets corresponding to left- and 
right-handed quarks and leptons. In this scenario tan /3, 
M H ±, and Msusy are free parameters. As far as CP 
phases are concerned, we adopt, without loss of gen- 
erality, the convention Arg(/i) = 0, while we assume a 
common phase for all the Af terms, &a = Arg(A f ) = 
Arg(-Ab) = Arg(A T ). As a consequence of this, we 
end-up with two free physical phases: <Pa and <P 3 = 
Arg(M 3 ). 

In addition to the parameters fixed by the CPX sce- 
nario, we need to fix the gaugino masses for our 
study. We take them as free parameters independently 
of M3 since, for them, we chose to relax the usual 
relations at the electro- weak scale: Mi/Mj — gf/g] 
with gi.j=gauge coupling constants, which originate 
from the assumption of gaugino-mass unification at 
the GUT scale. The neutralino x 1S defined as usual 
as the lowest- mass linear superposition of _B-ino B, 
VU-ino Vy( 3 \ and of the two Higgsino states H®, H®: 



X = a 1 B + a 2 W {3} + a 3 H° + a 4 H°. 



(6) 



In Ref. [4] it was proved that in a CP-conserving 
effective MSSM with \M X \ « |M 2 | light neutralinos 
of a mass as low as 7 GeV are allowed. Indeed, for 
I Mi I << I M2 1 the LEP constraints do not apply, and 
the lower bound on the neutralino mass is set by the 
cosmological bound. In the following we will assume 
vanishing phases for Mi and M2, and we will fix for 
definitencss M2=200 GeV (the phenomenology we are 
interested in is not sensitive to these parameters in a 
significant way). On the other hand, we will vary Mi, 
which is directly correlated to the lightest neutralino 
mass m x . 

In this work, we rely on CPsuperH [5] for the com- 
putation of mass spectra and couplings in the MSSM 
Higgs sector. 



3 Experimental constraints on the CPX 
scenario 

In the CPX scenario the lightest Higgs boson is mostly 
CP odd and its production at LEP is highly suppressed 
since |fl>HiVv| ^ 1 though it is kinematically accessi- 
ble. As a consequence of this, taking <P A = <l> 3 = 90° 
and Msusy = 0.5 TeV, the combined searches of the 
four LEP collaborations at y/s — 91 — 209 GeV re- 
ported the following allowed interval for a very light 
Hi [6], which we will focus on in our analysis: 
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M Hl < 10 GeV for 3 < tan/3 < 10. 



(7) 
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Fig. 1. The Thallium EDM d T i = dri X 10 24 tonfk in 
the CPX scenario with AfsusY = 0.5 TeV in the region 
Mh x 15 GeV and 3 < tan/3 < 10. The different shaded 
regions correspond to different ranges of |g!ti|, as shown: 
specifically, the narrow region consistent with the current 
thallium EDM constraint, \dri\ < L is denoted by black 
squares. In the blank unshaded region we have |g!ti| > 100. 
The region below the thick solid line is excluded by data on 
T(1S) decay. For comparison, the thin line shows an esti- 
mation of the same boundary obtained using the tree-level 
coupling taking O a \=l, i.e. Is^^J = tan/3. Also shown 
are the three contour lines of the reseated B(B S —* p/j) = 
(B s -> fifj) x 10 7 : B(B 3 -» nfj,) = 2 (solid), 20 (dotted), 
and 200 (dashed). 



We observe that in the scenario analyzed by the LEP 
collaborations one has |//|=2 TeV. For this large value 
of the neutralino is a very pure i3-ino configura- 
tion, with a Higgsino contamination a 3 ~0.02. As will 
be shown in Section[4l this has important consequences 
for the phenomenology of relic neutralinos, in partic- 
ular suppressing their annihilation cross section, and 
restricting the possibility of having a relic abundance 
in the allowed range only to the case of resonant an- 
nihilation. So, the exploration of different possibilities 
with lower values of \fj,\ could in principle be very rele- 
vant for relic neutralinos. However, this would require 
a re-analysis of LEP data which is beyond the scope 
of this paper. 

CP phases in the MSSM are significantly constrained 
by the EDM measurements. In particular, the EDM of 
the Thallium atom provides currently the most strin- 
gent constraint on the MSSM scenario of our inter- 
est. In Fig. [TJ we show the rescaled Thallium EDM 
dxi = ^Tl x 10 24 in units of earn in the M^-tan/3 
plane. Here, we consider only the contributions from 
the Higgs- mediated two-loop diagrams [7]. Different ranges 
of I (in I are shown explicitly by different shadings. In 
the blank unshaded region we have obtained |dn| > 
100. We also note that the Thallium EDM constraint 
can be evaded by assuming cancellations between the 
two-loop contributions and other contributions, such 
as those from first- and second-generation sfermions. 



Stefano Scopel 



The Lightest Higgs Boson and Relic Neutralino in the MSSM with CP Violation 



In this way the allowed region shown in Fig. Q] can 
be enlarged. The amount of cancellation can be di- 
rectly read-off from Fig. [T] For instance, in the region 
|dri| < 10 it would be less severe than 1 part in 10. 

In the region of Eq. (J7J) the bottomonium decay 
channel T(IS) — > jHi is kinematically accessible jS], 
and the experimental upper bounds on this process 
can be directly converted to a constraint in the plane 
tan /3-Mjii ■ The result is shown in Fig. [IJ where the 
thin (red) line corresponds to the limit obtained by ne- 
glecting finite-threshold corrections induced by gluino 
and chargino exchanges, and by setting O a \ = 1, while 
the thick solid line shows the same constraint when 
threshold corrections and the true value of O a i are 
used. From Fig. [T]one can see that, when the follow- 
ing constraints are combined: (i) the LEP constraint; 
(ii) Thallium EDM; (iii) the limit from bottomonium 
decay, the allowed parameter space is reduced to: 

7 GeV<M Hl <7.5 GeV and tan/3~3. (8) 

This region may be enlarged to 

7 GeV < M Hl < 10 GeV and 3 < tan/3 < 5 , (9) 

if we assume 10 %-level cancellation in the Thallium 
EDM. 

In light of the above discussion and for defimteness, 
from now on we will fix Mh 1 =7.5 GeV and tan/3=3 in 
our analysis. Taking into account the CPX parameter 
choice of Eq.© with $ A = $ 3 = 90° and M SUSY = 
0.5 TeV, this implies, in particular: Mjy± ~ 147 GeV, 
M H2 ~ 108 GeV, M H3 ^ 157 GeV. 

As will be discussed in the following sections, if 
the pseudoscalar Higgs boson mass is in the range (J5J) , 
a CPX light neutralino with m x < Mh ± /2 can be a 
viable DM candidate. Due to their very pure B-mo 
composition, and to the quite low value of tan (3, neu- 
tralinos in this mass range evade constraints coming 
from accelerators. For instance, in the CPX light neu- 
tralino mass range the present upper bound to the 
invisible width of the Z-boson implies \a 2 — a|| < a 
few percent, a constraint easily evaded in this case. 

A potentially dangerous constraint can come from 
the decay B s — > pp, since its dominant SUSY contri- 
bution scales as tan 6 (5 \p\ 2 /Mfj and may have a res- 
onance enhancement when H± is so light that Mh x ~ 
Mb s - In Fig. [IJ we show three contour lines of the 
rescaled B(B S — > pp) = (B s — ► pp.) x 10 7 : B(B S — ► 
pp) = 2 (solid), 20 (dotted), and 200 (dashed). For 
the parameters chosen by combining the results from 
LEP2 searches, Thallium EDM, and Bottomonium de- 
cay, Eq. we get: B(B S — > fifi)cpx — 6 x 10~ 7 
taking fg s — 0.23 GeV. This is three times larger 
than the present 95 % C.L. limit: B(B S — > pp) < 
2 x 10~ 7 . This can be easily made consistent with the 
present experimental constraint if some mild cancella- 
tion takes place. The "GIM operative point" mecha- 
nism discussed in Ref. [9] may be an example of such 
cancellation, when the squark mass matrices are flavour 
diagonal. In particular, we find that B{B S — > pp)cpx 



is consistent to the experimental upper bound by choos- 
ing 0.8 < p < 0.9, where p = mq/MsusY is the hierar- 
chy factor introduced in Ref. [9] , with mq the soft mass 
for squarks of the first two generations. 

For the discussion of other constraints, see [T]. 



4 The relic density 

Taking into account the latest data from the cosmic 
microwave data (CMB) combined with other observa- 
tions [TU] the 2-cr interval for the DM density of the 
Universe (normalized to the critical density) is: 



0.096 < Q m h z < 0.122. 



(10) 



where h is the Hubble parameter expressed in units 
of 100 km s" 1 Mpc -1 . In Eq.(fID| the upper bound 
on fi m h 2 establishes a strict upper limit for the abun- 
dance of any relic particle. In absence of some reso- 
nant effect, the natural scale of the annihilation cross 
section times velocity a ann v of CPX light neutrali- 
nos is far too small to keep the relic abundance be- 
low the upper bound of Eq. (|10jl (in particular they are 
very pure i3-inos and their mass is below the thresh- 
old for annihilation to bottom quarks, which is usu- 
ally the dominant channel of a ann v for light neutrali- 
ties [4]). However, when m x ~ M^/2 neutralinos 
annihilate through the resonant channel XX H\ — > 
standard particles, bringing the relic abundance down 
to acceptable values. In the Boltzmann approximation 
the thermal average of the resonant a ann v to the final 
state / can be obtained in a straightforward way from 
the following relation among interaction rates: 



< o- ann v > Ies j=< r(xx -> f) >= 



< r( xx -> H 1 )B(H 1 - /) >= 



n Hl r } 



■gl(gHi) 

K 2 {x Hl ) 



B 



(11) 



where brackets indicate thermal average, r x is the 
zero-temperature Hi annihilation amplitude to neu- 
tralinos and the thermal average of this quantity is ac- 
counted for by the ratio of modified Bessel functions 
of the first kind Ki and K 2 , Bf is the Hi branching 
ratio to final state /, n, = gim^K^Xi) / '{2-K 2 Xi) are 
the equilibrium densities with Xi = m^/T, T the tem- 
perature, and gi the corresponding internal degrees of 
freedom, g x = 2, gg^l. The factor of 1/2 in front of 
Eq. lfTTj) accounts for the identical initial states in the 
annihilation. 

The result of our calculation is shown in Fig. [2l 
where the neutralino relic abundance fl x h 2 is shown 
as a function of the mass m x . The asymmetric shape 
of the curve in Fig. [2] is due to the fact that thermal 
motion allows neutralinos with m x < Mh 1 /2 to reach 
the center-of-mass energy needed to create the reso- 
nance, while this is not possible for m x > M^/2. In 
the same figure, the two horizontal lines indicate the 
range of Eq. (p~0|) . 
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Fig. 2. Relic abundance 
mass m x for the CPX 
tan p=8, Msvsy=0.5 TeV and <P A = $3 = 90° 
horizontal lines indicate the interval of Eg. SKA) . 
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as a function of the neutralino 
scenario with Mh 1 =7.5 GeV, 
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that one can have < a ann v >o<C < o- ann v >. In fact, as 
already shown in Section [51 the thermal motion in the 
early Universe {x x ~ Xf ~ 20) allows neutralino res- 
onant annihilation when m x < Mh 1 /2. However, for 
the same neutralinos the contribution of the resonance 
to < a ann v >o can be negligible at present times, since 
their temperature in the halo of our Galaxy is of order 
x x ,o ~10~ 6 < ij. This implies that the annihilation 
cross section can be large enough in the early Uni- 
verse in order to provide the correct relic abundance, 
but not so large at present times as to drive indirect 
signals beyond observational limits. As a consequence 
of this, in [1] we have discussed in detail signals for in- 
direct Dark Matter searches and shown that they are 
compatible with the present experimental constraints, 
as long as m x < Mh x /2. On the other hand, part of the 
range m x > allowed by cosmology is excluded 

by antiproton fluxes. Moreover, in [T] prospects of de- 
tection in future DM searches have been discussed, 
showing that CPX neutralinos might indeed produce 
a detectable signal. Finally, as far as a very light Higgs 
boson Hi is concerned, we observe that the LHC might 
not be able to detect it, and a Super B factory could 
thus be needed for its observation fl"2l. 



In Fig. [2] the neutralino mass range allowed by cos- 
mology is: 3.15 GeV ~ m x ~ 3.83 GeV. Allowing for 
the variation of Mfjj within the range of Eqs. (|8I9[ ) , 
this range is enlarged to: 



2.93 GeV < m v < 5 GeV. 



(12) 



In this scenario the neutralino relic abundance can fall 
in the range of Eq. (fT0|) only with some level of tun- 
ing at the boundaries of the allowed mass range. For 
intermediate values of m x either the neutralino is a 
sub-dominant component of the DM, or some non- 
thermal mechanism for its cosmological density needs 
to be introduced. Of course all our considerations are 
valid if standard assumptions are made about the evo- 
lution of the early Universe (e.g. about the reheating 
temperature at the end of inflation, the energy budget 
driving Hubble expansion, entropy production, etc). 



5 Dark matter searches 

Neutralinos in the halo of our Galaxy can be searched 
for through direct and indirect methods. In particular, 
due to their very low mass and cross section, CPX light 
neutralinos are quite hard to detect through direct de- 
tection, although one proposal exists for such a hard 
task [TT] . As far as indirect searches of CPX light neu- 
tralinos are concerned, in our scenario the neutralino 
relic density Q x h 2 is driven below the observational 
limit by the resonant enhancement of the annihilation 
cross section < a ann v >. The same cross section calcu- 
lated at present times, < o~ ann v >o, enters into the cal- 
culation of the annihilation rate of neutralinos in our 
galaxy. This could produce observable signals, like 7's, 
i/'s or exotic components in Cosmic Rays (CR), like 
antiprotons, positrons, antideuterons. Note, however, 
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